Role of slip between a probe particle and a gel in microrheology 
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In the technique of microrheology, macroscopic rheological parameters as well as information 
about local structure are deduced from the behavior of microscopic probe particles under thermal 
or active forcing. Microrheology requires knowledge of the relation between macroscopic parameters 
and the force felt by a particle in response to displacements. We investigate this response function 
for a spherical particle using the two-fluid model, in which the gel is represented by a polymer 
network coupled to a surrounding solvent via a drag force. We obtain an analytic solution for 
the response function in the limit of small volume fraction of the polymer network, and neglecting 
inertial effects. We use no-slip boundary conditions for the solvent at the surface of the sphere. 
The boundary condition for the network at the surface of the sphere is a kinetic friction law, for 
which the tangential stress of the network is proportional to relative velocity of the network and the 
sphere. This boundary condition encompasses both no-slip and frictionless boundary conditions as 
limits. Far from the sphere there is no relative motion between the solvent and network due to the 
coupling between them. However, the different boundary conditions on the solvent and network tend 
to produce different far-field motions. We show that the far field motion and the force on the sphere 
are controlled by the solvent boundary conditions at high frequency and by the network boundary 
conditions at low frequency. At low frequencies compression of the network can also affect the force 
on the sphere. We find the crossover frequencies at which the effects of sliding of the sphere past 
the polymer network and compression of the gel become important. The effects of sliding alone can 
lead to an underestimation of moduli by up to 33%, while the effects of compression alone can lead 
to an underestimation of moduli by up to 20%, and the effects of sliding and compression combined 
can lead to an underestimation of moduli by up to 43%. 

PACS numbers: 83.10.-y, 82.70.Gg, 83.50.Lh, 83.60.Bc 



Single particle microrheology [i|, 0] has been a useful 
tool for measuring the material characteristics in situ- 
ations where traditional rheometers are difficult to use. 
For example, it has been particularly useful in obtaining 
rheological measurements when large sample sizes are dif- 
ficult to obtain, and when removing samples from their 
natural environment may be undesirable, such as in living 
cells. Traditional rheological measurements typically ob- 
tain the frequency-dependent macroscopic shear modulus 
G{uj), which relates shear stresses to homogeneous shear 
deformations. To interpret the results of single particle 
microrheological measurements, one must know what the 
response of a particle embedded in the material of ques- 
tion will be to driving forces (or equivalently, the force 
exerted on the particle in response to particle displace- 
ments). In practice, the response of such a particle has 
been assumed to relate to the macroscopic shear modulus 



-67raRe{G'(w)(5rexp(— loji)} . 



(1) 



when the particle with radius a oscillates at frequency 
w and displacement Re{(5r exp(— iwt)}. The modulus 
G = G' — iG" is in general complex. In microrheological 
experiments, the use of this force response in conjunction 
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with the fluctuation-dissipation theorem has been called 
the "generalized Stokes-Einstein relation" |l| . 

Eq. [1] holds true when the material may treated using 
continuum mechanics as a single incompressible phase 
with a complex shear modulus and no-slip boundary con- 
ditions between the particle and the material 2, 3]. The 
generalized Stokes-Einstein relation has been validated 
experimentally in the test case of a solution of wormlike 
micelles [J] , but in many cases there are a number of is- 
sues that can complicate microrheological measurements, 
including compressional effects [2, 5], local depletion of 
the polymer [y, |7|, modification of local properties via 
surface chemistry of the particles @, Q i and violation of 
no-slip boundary conditions [lOj . 

Although the technique of 2-particle microrheology p, 
ill| ameliorates many of these problems, it is much more 
technically demanding than single-particle microrheol- 
ogy. Furthermore, fully understanding how all the above 
effects can impact 1-particle microrheology allows us to 
use it to understand materials and processes in which lo- 
cal heterogeneity plays an important role. In this paper 
we demonstrate how to quantify the effects arising from 
violation of the no-slip boundary condition for both an in- 
compressible and compressible gel. Previous studies have 
treated the response of the medium to aparticle in the 
case of no-slip boundary conditions [1, y, H, [l2l and in 
the presence of compression [5| . In this work we include 
the effects of sliding between the sphere and the medium 
in a two-fluid model for a gel with phases representing a 
viscoelastic polymer network and fluid solvent. Through- 



out this paper, we work in the dilute network hmit, which 
is appropriate for many microrheological studies, such as 
those of actin networks @ and DNA solutions [3| . Slid- 
ing is implemented using a kinetic friction law — the shear 
stress of the network is proportional to the relative ve- 
locity between the network and the surface of the sphere, 
with proportionality constant S. No-slip (S = oo) and 
frictionless sliding (S = 0) boundary conditions can be 
obtained as limiting cases. 

Physically, we can see the effects of slip by comparing 
two different situations: a solid sphere of radius a moving 
with velocity v through a liquid of viscosity 77, and a clean 
bubble of radius a moving through the same liquid. In 
the first case, there is no slip at the surface of the sphere, 
and the drag force is f = —Gnri'va. In the second case, 
the drag on the bubble is smaller, f — —ATTrj-va. In the 
far field, the velocity field for the fluid flow for the bubble 
is also reduced by a factor of 2/3 relative to flow for the 
solid sphere |13| . 

Our use of the two-fluid model for a gel allows us 
to include the effects of sliding via a boundary condi- 
tion between the sphere and polymer network, even as 
the solvent retains no-slip boundary conditions. Because 
the fluid and network are coupled by drag, far from the 
sphere there is no relative motion between the network 
and the solvent. However, because the fluid and network 
have different boundary conditions, the far-field motion 
can have the character of that driven by no-slip boundary 
conditions, that driven by frictionless sliding boundary 
conditions, or somewhere in between. We find that the 
far-field motion is controlled by frequency. In the high 
frequency limit, the far- field solution has the properties of 
a deformation driven by no-slip boundary conditions. In 
this case, the solvent flow is the same as the flow around 
a sphere in a simple viscous fluid and no-slip boundary 
conditions, and the network is dragged along by the sol- 
vent. In the low frequency limit, the far- field solution 
has the properties of a deformation driven by frictionless 
sliding boundary conditions. The network displacement 
is the same as the displacement around a sphere in a sim- 
ple elastic solid and frictionless boundary conditions, and 
the drag prevents the solvent from moving faster than the 
network. 

Similarly, the force felt by the sphere interpolates be- 
tween the limits provided by no-slip and frictionless slid- 
ing boundary conditions. The effects of sliding and com- 
pression can be described using three crossover frequen- 
cies. In the following, a is the radius of the probe particle, 
I is the mesh size of the polymer network, 77 is the viscos- 
ity of the solvent, fi and A are effective Lame coefficients 
of the network (Eq. [2]), and S is the friction coefficient 
between the probe particle surface and the polymer net- 
work (Eq. [Ti)) . 

1. Associated with sliding, ujs ~ \^\l/{r]a) for a > I, 
and LPs ~ \^\l'^ / iv^^) for a < I. 

2. Associated with sliding friction, LUf w |/i|/(Sa). 

3. Associated with compression, Wc ~ \2fi+X\P/{ria^). 



The properties of the force can be summarized as 

(a) At frequencies larger than lj^, w/, and tOc, the re- 
sponse force obeys Eq. [U and moduh deduced us- 
ing the generalized Stokes-Einstein relation match 
macroscopic measurements. 

(b) At frequencies smaller than lUs and ujf, but larger 
than ujc, effects from sliding are important, the re- 
sponse force is reduced, and moduli deduced using 
the generalized Stokes-Einstein relation underesti- 
mate by up to 33%. 

(c) At frequencies smaller than Uc, but larger than Ug 
or ujf, effects from compression are important, the 
response force is reduced, and moduli deduced us- 
ing the generalized Stokes-Einstein relation under- 
estimate by up to 20%. 

(d) At frequencies smaller than Ug, w/, and Wc, effects 
from both sliding and compression are important, 
the response force is reduced, and moduli deduced 
using the generalized Stokes-Einstein relation un- 
derestimate by up to 43%. 

(e) In real systems, fi and A may be frequency de- 
pendent, and therefore so are Ws,/,c(w). In light 
of the above, if at any point in an experiment 
Lo < u!sj,c{'-^), caution should be exercised in apply- 
ing the generalized Stokes-Einstein relation. Our 
results provide a way to analyze microrheological 
data in this more complicated situation. 

The paper is organized as follows. In Section[T]we intro- 
duce the two- fluid model for a gel and discuss how sliding 
boundary conditions can lead to qualitatively new behav- 
iors. We discuss these behaviors in the simple geometry 
of a gel between two oscillating plates. In Section [III we 
describe the analytic solution for the force exerted on an 
oscillating sphere in a gel, with some details relegated to 
Appendix]^ Then we describe the properties of the force 
in Section imi All the features associated with sliding are 
present in the case of an incompressible network, which 
we discuss in the main body of the paper. In the dis- 
cussion we describe the implications for microrheological 
experiments, providing examples of how our results can 
be used to help interpret microrheological data. Finally, 
in Appendix |B] we describe the properties of the force in 
the presence of both sliding and a compressible network 
for completeness. 



I. THE IMPORTANCE OF BOUNDARY 

CONDITIONS IN THE TWO-FLUID MODEL 

FOR A GEL 

We describe the gel as a polymer network with Lame 
coefficients fj, and A, interspersed with fluid with viscosity 
77 (Fig. [1} . In general, we assume that the network shear 
modulus is complex, and write n = /ui— z/i2- The polymer 




FIG. 1; A probe particle modeled as a sphere of radius a 
oscillating with frequency li; in a gel modeled by two phases 
representing a viscous solvent and a viscoelastic polymer net- 
work with mesh size I. 



and fluid are coupled to each other by a friction coefficient 
r [^, HJ, \1M, ■ We ignore any inertial contributions to the 
equations of motion: 



/^V2u+ (A + /i)V(V-u) = r(u-v) 

77VV - Vp = -r (ti - v) 

V • [(1 - ip)v + ipu] = 0. 



(2) 
(3) 
(4) 



Here u is the displacement field of the polymer network 
and V is the velocity field of the solvent, (p is the volume 
fraction of the polymer network, and Eq. |4] expresses the 
overall mass conservation of the gel. In the limit of small 
volume fraction, we can approximate Eq.[4]as the incom- 
pressibility of the solvent, V • v = 0. In Eqs. [2] and [3 
the dot denotes a material time derivative; throughout 
this paper we replace this with a partial time derivative 
since we only consider small displacements and work with 
these equations only to linear order. To deduce the forces 
exerted by the gel, we use the stress tensors 



rPoly 



/i [Vu + (Vu)'^] + I AV • u 






(5) 

(6) 

(7) 



It is important to note that the moduli are effective mod- 
uli for the network in the presence of the solvent includ- 
ing, for example, osmotic effects. In addition, fi and A 
have implicit dependence on volume fraction (f) since the 
network stiffness depends on the network density. For a 
typical microrheological experiment, we are interested in 
a dilute network, with mesh size much larger than the 
network filament diameter. For such a dilute network, 
r « i^/P, where I is the mesh size of the network Q. We 
note that Eqs.[2]and[3]are not restricted to the dilute net- 
work limit; for example, in a dense network, these equa- 
tions reduce to the poroclastic equations if x*^""* is neg- 
ligible compared to the other stresses in the system [l^ . 



This model has previously been used by Levine and 
Lubensky [5| to confirm that even with compressional 
effects, the force on a sphere obeys Eq. [1] for a range 
of frequencies. Above this range of frequencies inertial 
effects become important, while below this range of fre- 
quencies, the compressibility of the material can become 
important. So far. the effects of compressibility have not 
been observed |17| . 

In this work, we also use the two-fluid model, so we can 
incorporate the compressibility effects described above. 
We ignore inertial effects, which typically are not im- 
portant at frequencies used in microrheological measure- 
ments. However, while previous results were obtained for 
no-slip boundary conditions, we obtain exact analytical 
results using boundary conditions which allow the sphere 
to slide with respect to the polymer network. Our bound- 
ary conditions encompass both the no-slip and friction- 
less sliding limits with respect to the polymer network. 

Although Norris [l^l and Oestreicher ^3] have previ- 
ously obtained analytic solutions for the response func- 
tion of a sphere in a material with a single complex shear 
modulus for both no-slip and frictionless sliding bound- 
ary conditions, a crucial feature of the two-fluid model is 
that the boundary conditions for the fluid and the net- 
work can be different. This feature opens up the pos- 
sibility that the far-field behavior can reflect either the 
fluid or the network boundary conditions, depending on 
the frequency of motion. 

A simple demonstration of this frequency-dependent 
behavior can be seen in a gel between two oscillating 
plates located at y — ±L/2 (Fig. [2^). The plates os- 
cillate with amplitude ±b cos ujt along the x direction, 
with no-slip boundary conditions between the plates and 
the solvent, Vx{y — ±L/2) = ^bujsinujt. First consider 
no-slip boundary conditions between the plates and the 
network, Ux{y = ±L/2) = ±bujshiujt ^]. In this case 
there is no relative motion between the network and the 
solvent, and both undergo simple shear motion for all 
frequencies. Since the strain is homogeneous, the stress 
exerted on upper plate is by definition determined by 
the macroscopic shear modulus G. Writing this stress as 
o'xy = Re[CTexp(— zwi)] leads to 



G 



a 



b/{2L) 



Ml -«(M2 + iiuj). 



(8) 



When the network can slip past the plate, the situation 
becomes more complicated, since the strain need not be 
homogeneous. For example, suppose there is frictionless 
sliding between the plates and the polymer network, so 
that the tangential stress on the network vanishes at the 
plates, dyUx{y = ±L/2) = 0. For simplicity assume /i is 
real (/i2 = 0) . Only the x-components of the velocity and 
displacement fields are nonzero, and the solutions are 



-Re 



-Re 



■2^bu^ycosh!^ 

kLrju) 



- sinh ky J 



— 2/ cosn — + smn Ky I 

77 Tr exp (—iLut) > , 

cosh^ + 2sinh^ J 

2ibuj (^y cosh !^-!f sinh ky) 

^ exp {—ibjt) 






cosh ^ + 2 sinh ^ 



(9) 
(10) 



r 



In these equations we introduce the complex quantity 
k = •\/r(l/?7 — iuj/ fi) which impUes an associated length 
scale, 



d = 



{v/txTaFWa^} 



Re 



(11) 



The penetration depth d determines the thickness of the 
layer next to each plate where the solvent moves rela- 
tive to the network. Friction suppresses relative motion 
beyond the penetration depth. 

The velocity and displacement fields are plotted in 
Fig-Hh- To see the range of behavior that is possible, it is 
useful to consider limiting cases. First, suppose there is 
no coupling between the network and the solvent, F = 0. 
Then the network has zero displacement and the solvent 
has homogeneous oscillatory strain (similar to solid black 
lines. Fig. [5]). Now suppose F 7^ 0. At sufficiently high 
frequency, the penetration depth d ~ •\/2/i/(Fa;) is small 
compared to L, and there is no relative motion between 
the network and the solvent in most of the gel. Deep 
within the gel, the solvent has a constant strain rate and 
the network has a constant strain. The penetration depth 
d determines the size of the boundary layers where the 
displacement changes from uniform strain to no strain 
to satisfy the condition of zero tangential stress at each 
plate. On the other hand, when the frequency is suffi- 
ciently low, d is approximately the mesh size I and we 
may still have L ^ d if F is big enough. In this case, 
again there is no relative motion between the solvent and 
the network deep in the gel, but now the deformation of 
the network is small. The penetration depth determines 
the thickness of the boundary layers at each plate where 
the velocity of the solvent rises from the small velocity of 
the network to the velocity of the plate. 

To determine the crossover frequency that marks the 
transition between the high- and low-frequency behav- 
ior, it is convenient to plot the strain of the network at 
the midplane of the gel, dyUxiO), as a function of dimen- 
sionless frequency, tjuj / ^l (Fig. [2};) . At high frequency, 
we see that dyUx{Q) ~ h/2L - the network is dragged 
along with the solvent and the motion in the interior of 
the gel is simple shear. At low frequency, ^^^^^(O) ~ 
- the deformation of the network and velocity of the 
solvent in the interior of the gel are small. Figure [5J: 
shows the crossover for various values of the coupling F. 
The scaling of the crossover frequency cjpiato with F can 



be deduced by examining the spacing of the curves in 
Fig[2t. For F > rj/L^, the crossover frequency is approx- 
imately Wpiatc ~ l^/{I^V^): with power law F~^'^; while 
for F < Tj/Lp' , the crossover frequency is approximately 
Opiate ~ /^/(^^r); with power law F~^. Later we will see 
similar behavior of the crossover frequency as a function 
of F for the spherical geometry. 

When F > rj/L^ , we observe from Fig. [2j: that the 
crossover frequency satisfies r;cijpiatc//^ ^ 1, and thus con- 
clude that d ~ \/v/^ loi' '^ near cjpiato- Let us further 
suppose that d <^ L. In this case we can understand 
the crossover frequency by examining force balance on 
the gel. As emphasized earlier, the solvent and network 
move together except in the boundary layers near the 
plates. The plates exert no direct force on the network, 
so any shear on the network results from the motion of 
the solvent dragging the network along in the bound- 
ary layer, and any stress from network shear in the in- 
terior is ultimately supported by the traction of the sol- 
vent at the plate. The forces can be estimated by ap- 
proximating Vx and Ux with piecewise linear functions, 
as in Fig. [31 In the limit of zero Reynold's number, 
there is force balance on any given layer of gel. Con- 
sider a layer bounded by the plate on one side, with 
width d. At the plate, the traction is of magnitude 
77(^(6 — Ue)/d, where Ue is the displacement at distance 
d from the plate. At the other side of the layer, the 
solvent and network stresses add up to approximately 
{riLu + fi)ue/{L/2 — d). Equating these two stresses deter- 
mines Ue- Using dyUx{0) « Ue/{L/2 — d), we find that the 
normalized strain dxu{0)/[b/{2L)] « {l/[l + 2fj.d/{r]Luj)]. 
For to > ^id/{r]L) ^ fj,/{Ly/Trj) the normalized strain 
is approximately one, meaning that u ~ b, while for 
Lu < ^./{Ly/Tri) the normalized strain decreases rapidly. 
Therefore, our force balance argument yields the same 
crossover frequency we identified from Fig. [2j;. Physi- 
cally, the traction due to the network shear is supported 
by the drag force between the network and solvent in 
the boundary layer of relative motion. The higher the 
frequency, the faster the relative motion, the more drag 
force, and the greater the network shear will be in the 
interior of the gel. 

This crossover behavior also shows up in the stress ex- 
erted by the gel on the plates. Suppose the force on the 
plates is interpreted as a measurement of the macroscopic 
shear modulus. Again expressing the stress exerted on 
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FIG. 2: (color online) a) Dimensionless amplitude of oscilla- 
tions as a function of dimensionless height for solvent velocity 
(thin lines) and polymer network displacement (thick lines) 
in a gel between two oscillating plates. The polymer network 
has frictionless boundary conditions at the plate surface. The 
dashed lines have VL^ /rj = 10^'^ and L <^ d. Far away from 
the plate the network and solvent move together. The solid 
line has FL^ /rj = 10"^'^, rJ^/^J■ = 0.01, and d :$> L. In this 
case the network and solvent need not move together as a unit 
far from the plate, b) Semilogarithmic plot of the dimension- 
less slope of the amplitude of oscillations of the network dis- 
placement at J/ = as a function of dimensionless frequency. 
The shape and spacing of the plots indicates that the crossover 
frequency tjpiato behaves as r~^'^ for VL^ /r) > 1, and as V~^ 
for TL'^/ri < 1. 



the upper plate as cFxy — Re[(Texp(— zcji)], we obtain 

Re a 



^P'"^"^ b/{2L) 






Im a 

b/Wr 



(12) 
(13) 



Due to the slip and relative motion, the gel does not 
undergo a homogeneous shear deformation, and GpcS 
need not be equal to G. 




FIG. 3: (color online) Simplified profile of dimensionless oscil- 
lation amplitudes for solvent velocity (blue dashed) and net- 
work displacement (red solid) for estimating force balance. 
The thin black line is the profile for a single fiuid between the 
plates. 



At high frequencies, the network is pulled by drag 
to move like the uncoupled solvent, and G' g ~ ^i 
(Fig. [3^), showing the response of an elastic solid be- 
ing sheared between the plates. Below the crossover fre- 
quency, G' gjj decreases rapidly, taking the character of 
a network with frictionless sliding boundary conditions, 
which does not exert stress on the plates. Similarly, at 
high frequencies, G'^^^ff ~ ijuj (Fig. [4)d), showing the vis- 
cous response of a simple fluid with viscosity rj sheared 
between the two plates. In Fig. [JJa, the quantity plotted 
is G'' ^g/rjLu, which is proportional to an apparent viscos- 
ity of the medium. Below the crossover frequency and for 
L ^ d, the apparent viscosity is enhanced. The profile 
of Fig. [3] explains this enhancement: at low frequencies. 
Me is nearly zero, and the solvent shear is approximately 
buj/d. Therefore the viscous stress on the plate is ap- 
proximately rjbuj /d ~ F""^'^. In the rest of this paper we 
describe analogous physics which arises in the spherical 
geometry, and discuss its ramifications for microrheolog- 
ical experiments. 



II. SOLUTION FOR SPHERE MOVING IN GEL 

Consider a sphere of radius a surrounded by a gel. We 
solve for the flows and displacements in the frame of the 
sphere, using spherical coordinates (r, 9, 0). We consider 
a sphere oscillating along the z direction [6 = 0) with 
amplitude Bcos{Lut). 

In the frame of the sphere, the boundary con- 
ditions for the fluid and solid at r — oo 
are v(oo,6', (/)) = — (cos^, — sin^, 0)i3a;sin(a;t) and 
u((X),6', (/)) — {cos9, — sm0,O)Bcos{ujt). The boundary 
conditions at the surface of the sphere are 



v{a,9,( 
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(14) 
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FIG. 4: (color online) a) Semilogarithmic plot of dimension- 
less Gp ;,fi determined by the stress exerted by the gel on the 
plates, as a function of dimensionless frequency, b) Loga- 
rithmic plot of G'p,,^/{rjuj) as a function of dimensionless fre- 
quency. Gpcff is scaled by the frequency so the plotted quan- 
tity is proportional to an apparent viscosity. 



Here the dot denotes a time derivative, and crP°'^™'^' is the 
stress tensor of the polymer network. The last condition 
allows the polymer network to slide past the surface of the 
sphere in the tangential direction, with a stress exerted 
proportional to the relative velocity. The magnitude of 
the force is controlled by the friction coefhcient S. For 
S ^ cxD we obtain no-slip boundary conditions, while for 
S = we obtain frictionless sliding. 

In dimensionless form, Eqs. HJHand Eq. [HI are 



— curl^u V(V • u) — — 7(u — v) 

if p 

curl^v-fVP ^ 7(u-v) 



(15) 
(16) 



v(cx3,6', 0) = — (cos0, — sin0,O)&sini 

u((X),0, 0) = (cos 6*, — sin 6*, 0)6 cost 

v(a,6'(^) = 

Ur{a,9,4>) = 

aP°"'"^'='-(a,0,0) = C^«(a,0,0), (17) 

where we measure lengths in units of a, time in units 
of 1/w, and pressures and stresses in units of rjoj. The 
dimensionless parameters are R = tjlu/ p., p = rjoj I {2p + 



A), 7 = To? /rj, b = B/a, and ^ = Ea/rj. As mentioned 
before, F « rj/l^ 0], where I is the mesh size of the 
polymer network, so that 7 w a^/P. For convenience 
we keep the same notation for dimensionless fields as we 
used for dimensional fields. 

Our method is an extension of the classic solutio n by 
Stokes of a sphere moving through a viscous fluid [l3|. 
The velocity field is divergence-free and can be expressed 
in terms of a single stream function $: 



curl 



^0 



(18) 



On the other hand, the polymer is compressible, so the 
displacement field has both a stream function for the 
divergence- free part ($) , and a potential function for the 
curl- free part (x): 



curl 



$0 



Vx- 



(19) 



Taking the curl and divergence of the above equations, 
leads to 



-v^x = -tV^x 
p 



(20) 



-7CurW,^-,^)(21) 
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The solution for x is written in the appendix, and in- 
volves three undetermined coefficients. 

Diagonalizing the equations for $ and ^ yields 




curl 



where 



a \ <p 
P J r sin ( 



— curl 
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P I r sm ( 



M 



(23) 



(24) 



and M is a two-by-two matrix. The appendix contains 
expressions for the matrix M and eigenvalues Aq, Xp 
which diagonalize the equations. The appendix also con- 
tains the solutions for a and (3\ each involves three un- 
determined coefficients. 

The nine undetermined coefhcients are determined by 
imposing the boundary conditions at infinity and the 
sphere surface. In addition, one must demand that Eq. 
llSl is satisfied by the combination of stream and potential 
functions in Eq. [121 The conditions on the coefficients 
are displayed in the appendix and determine the velocity 
and displacement fields. 

From the velocity and the displacement fields the 
(nondimensional) stresses and pressure of the polymer 



and fluid can be calculated 
1 
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.fluid 



^[Vu+(Vun+I Q-i)V.u(25) 



The total stress tensor is 



total poly 



„fluid 



pi. 



(28) 



Vv + (Vv)"^ 
dr ['y{ur - 



W^Vr] 



(26) Integrating the stress tensor over the surface of the sphere 
yields the total force on the sphere, which is (in dimen- 

(27) sionless form) 
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In these expressions, whenever a square root appears, 
the root with positive real part should be chosen. This 
expression is one of the main results of this paper. Al- 
though it is in analytic form, it is still unwieldy. In the 
next section we summarize its properties. 



III. PROPERTIES OF RESPONSE FORCE 

The dimensional response force is Eq. [30] multiplied 
by a^rjuj. If interpreted using the generalized Stokes- 
Einstein relation, this yields effective moduli 



G„ff = — - — Re 
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If Eq. [T] holds, we would expect G'^g — fii, and G"ff = 
/i2 + tor). In the following we will describe the response 
force in terms of the effective moduli. The descriptions in 
the following sections are obtained by analyzing the form 
of Eq. [3D1 To give clarity to the discussion, in all the plots 
we assume that rj and p are frequency-independent. 

To understand the properties of the response force 
when the particle can slide past the network, first we look 
at the limit of incompressible network, A ^ oo, which 
displays all of the effects due to sliding boundary con- 
ditions. Network compressibility (finite A) complicates 
the situation by adding additional features to the effec- 
tive moduli. We leave a detailed description of the more 
complicated scenario with effects from both sliding and 
compressibility to the appendix, but summarize the re- 
sults in the discussion. 



A. Frictionless limit with incompressible network 

To illustrate the effects of sliding when the network 
is incompressible, we concentrate on the case of friction- 
less sliding between the sphere and the polymer network, 
^ = 0. When both the network and solvent are incom- 
pressible and have no-slip boundary conditions, there is 
no relative motion between the network and the solvent, 
and Eq. [1] holds [3]- In contrast, when the polymer 
network has frictionless boundary conditions at the sur- 
face of the sphere, the network and solvent move relative 
to each other near the sphere. This relative motion be- 
tween the network and solvent leads to drag forces which 
tend to diminish the relative motion, so that far from 
the sphere the network and solvent move together. The 
far field solution can be either that of an incompressible 
material driven by no-slip boundary conditions, that of 
an incompressible material driven by frictionless bound- 
ary conditions, or somewhere in between. As in the ex- 
ample of a gel between two oscillating plates, at high 
frequencies, the solution driven by the velocities (the sol- 
vent, with no-slip boundary conditions) wins out, while 
at low frequencies the solution driven by displacements 
(the polymer network, with frictionless boundary condi- 
tions) wins out. 

In the limit of frictionless sliding (^ ^ 0) and inconi- 
pressibility {p -^ 0) Eg . [501 reduces to 
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The low frequency limit is obtained when the first term 
in the brackets in the denominator dominates over the 
second term, and the high frequency limit is obtained 
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FIG. 5: (color online) Logarithmic plot of dimensionless G^ff 
as a function of dimensionless frequency with frictionless 
boundary conditions between sphere and an incompressible 
polymer network. Eq.[T] would predict G'^g — fii. The shape 
and spacing of the plots indicates that there is a crossover 
frequency ujs « |p|/(^a/7) for 7 > 1, and uJs « |/i|/('77) for 
7 < 1. For ij < UJs, effects from sliding cause G'^g to underes- 
timate /ii by up to 33%. 



in the opposite case. Therefore the crossover frequency 
LUs satisfies \Ry/^/^/l - iR\ w 1. For 7 sa a^/P > 1 
(sphere larger than mesh size), Us ~ ImI/('7\/7)i while 
for 7 < 1 (mesh size larger than sphere), tOg ~ iMl/l'yT)- 
The expression for oug in the introduction is written in 
terms of variables which can be easily controlled experi- 
mentally using the approximation 7 « o? /P, and hence 
applies when T « rj/P. Note that in these two limits, 
the form of the crossover is the same as in the case of 
the plates but with the length scale L replaced by a; in 
particular the dependence on F is the same as in the case 
of the oscillating plates, with power law T~^/'^ and T~^, 
respectively. 

For frequencies a; ^ Ws , the modulus GcS obeys Eq. [U 
so that GJjff = /ii and G"g = — (?yw -I- ^2)- For uj <C Wg, 
Goff behaves as 



Goff(w < Us 
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For the lowest frequencies, the real modulus G^g = 
2/ii/3, 67% of the value in Eq. [H The behavior of G;,^ 
as a function of frequency is plotted in Fig. [51 

Now we turn to the imaginary modulus G"jf . Eq. [1] 
predicts that G"g = /i2 + tjoj. In the example of the 
plates, we assumed /12 = 0, and found an enhancement 
in the effective viscosity at low frequencies. In the case of 
a sphere there is a similar enhancement to the effective 
viscosity, but sliding also affects the estimation of ^2- 

First consider the case where the network shear mod- 
ulus is purely real {1x2 =0). In this case, all dissipation 
arises from the viscous solvent and friction between the 
solvent and the polymer network. The behavior of G"jj 
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FIG. 6: (color online) Logarithmic plots of dimensionless 
G"ff as a function of dimensionless frequency for frictionless 
boundary conditions between the sphere and an incompress- 
ible polymer network, a) G'^^/{r]U)) as a function of frequency 
for /.12 = 0. This contribution to G"ff from solvent viscos- 
ity and network-solvent drag dominates for small ^2 or high 
frequencies. Inset: G"fi rescaled by its low frequency magni- 
tude (see text) to demonstrate that the crossover frequency 
is UJs- b) G"ff as a function of frequency. Eq.[T] would predict 
G"ff = fi2+'riuj. For uj < uJs, effects from sliding cause G"fi to 
underestimate fi2 by up to 33%. The crossover behavior is the 
same as that of the real part (G'^g) except that for u > 112/ri 
the solvent viscosity (rjui) dominates. 



with /i2 = is shown in Fig. ^i, G'^g/{riuj), an appar- 
ent viscosity. When uj '^ uj^, the apparent viscosity is 
ry. For frequencies uj <ti Ws, the apparent viscosity of the 
"solvent" contribution tends towards a plateau with en- 
hanced viscosity « 77(1 + y/j/9). Physically, as in the 
case of the plates, the solvent moves more than the net- 
work in a layer near the sphere. A region of enhanced 
solvent shear and viscous drag is set up as the solvent 
reduces its motion relative to the network, until far away 
the solvent takes the same motion as the network. As in 
the case of the plates, the enhancement depends on F as 
r^^/^. In the inset to Fig. [5^ we plot G"g rescaled by its 
low frequency value, 1 -I- (G"jj — riuj)/{r]ujy/^/9), to show 
that the crossover value is Us- 

Now consider the case when the imaginary part of the 
network modulus {(12) is nonzero. In addition to the con- 
tribution to G"g described in the previous paragraph. 



there are also contributions from the network modulus 
/i2- For w ^ Ws, the loading contribution is /i2, the re- 
sult of Eq. [TJ When lo ^ lUs, the leading contribution 
is 2/i2/3, 66% of the high-frequency contribution. For 
a constant /i2, at high frequencies the contribution pro- 
portional to "qLO dominates, while at low frequencies the 
contribution from /i2 dominates. In total, G'l^ behaves 
as shown in Fig. ^p. 



B. Incompressible network with intermediate 
friction 



For any finite value of ^, there will be a crossover from 
behavior similar to ^ = at low frequencies to behavior 
similar to ^ ^ oo at high frequencies, as shown in Fig. [T] 
The crossover is set by the condition \R\£^ = 1 and occurs 
at LOf fa I ^ I ■ Because the no-slip limit of G'^g is /ii , 
G'gff appears to crossover from a low frequency limit of 
2/ii/3 to the high frequency limit iii ai lu — min{a;y , ujs}- 
Similar behavior is seen in the imaginary modulus G"g. 



IV. DISCUSSION 

The main aim of this paper is to provide an under- 
standing of how sliding between a probe particle and 
polymer network can affect the response function in a gel. 
Motivated by experimental data, Starrs and Bartlett [10| 
have previously suggested the possiblity of frictionless 
sliding boundary conditions being relevant in microrhe- 
ological experiments. Our solution shows how this may 
come about, despite the fact that the liquid solvent al- 
ways couples to probe particles with no-slip boundary 
conditions in continuum fluid mechanics. We find that, 
due to the coupling between the solvent and polymer net- 
work, in far-field, the solvent can move as if it is driven 
by sliding boundary conditions. Specifically, when the 
fluid and polymer network have different boundary con- 
ditions, at high frequencies this coupling tends to lock 
the network into the motion of the fluid, while at low 
frequencies the coupling locks the fluid into the motion 
of the network. 

Our analytical solution for the force on a probe parti- 
cle in response to oscillatory motion quantifies the effect 
of sliding. We have described the effects of sliding on 
the force felt by probe particles in terms of an effective 
modulus Gcff . For high enough frequencies, GcS matches 
the amcroscopic shear modulus G, while for low frequen- 
cies Gcff may underestimate G. We have identified two 
crossover frequencies, one associated with sliding (lUs), 
and one associated with friction (w/). For any finite 
amount of surface friction, when frequency is reduced 
below both tUg and LUf {uj <C minjws, wy}), Gofr tends to 
underestimate moduli by 33%. 

The description in the previous paragraph applies in 
the case of an incompressible network, as treated in the 
body of this paper. However, in real materials the com- 
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FIG. 7: (color online) Logarithmic plot of dimensionless a) 
G^ff, b) G"ff as a function of dimensionless frequency for an 
incompressible network and for boundary conditions with in- 
termediate friction between sphere and polymer network. In 
these plots /ii = /i2 and 7 = 10*. The behavior interpolates 
between that of no-slip boundary conditions at high frequen- 
cies and frictionless boundary conditions at low frequencies, 
with crossover frequency ujf « \fi\/{Sa) . For frequencies be- 
low both u!s and ujf the moduli can be underestimated by up 
to 33%. 



pressibility of the network can lead to additional reduc- 
tions in Gcff. These compressional effects complicate the 
behavior of Gcff as a function of frequency, but may be 
relevant in experimental situations, so we have described 
them in detail in Appendix [Bl Intuitively, a compress- 
ible network is softer than an incompressible network and 
should exert a smaller force on an oscillating sphere. Our 
analytic solution identifies a third crossover frequency as- 
sociated with compressive effects (wc). This is the same 
crossover frequency identified earlier by other investiga- 
tors 0, 0- ^^ terms of the effective modulus Gcff, the 
effects of compression and sliding together are as follows: 
For frequencies much greater than ujs, ujf, and to a GcS 
matches the macroscopic modulus G. For frequencies 
such that Taiii{u!s,i-Uf} <^ uj <^ ujc, compressional effects 
reduce Gcff by up to 20%. For frequencies such that 
cjc ^ w ^ inm{LUs,^f}, sliding effects reduce Gcff by 
up to 33%. Finally, for the lowest frequencies, much less 
than LOs, LOf, and tOc, both compressional and sliding ef- 
fects reduce Geff by up to 43% relative to G. 
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For any specific experiment, our results describe how 
the appropriate effects of sUding on single-particle mea- 
surements should be accounted for. For example, if 
7 Ki a? /P > 1, which occurs for particles larger than 
the mesh size, then a good estimate for w^ is given by 
rju)s/\ij\ ~ l/-\/7 ~ l/o-- Therefore iUs can be estimated 
in experimental measurements, allowing determination of 
whether it is possible for corrections to Eq.[T]from sliding 
effects to be important. In general, /i will be frequency 
dependent, and therefore so will a;g(a;). A simple test 
which follows from our analysis is that if throughout the 
range of experimental frequencies to ^ Ws(w), then one 
can determine that sliding will not affect the effective 
moduli. In the following two paragraphs, we apply this 
test to two sets of experimental results reported in the 
literature. 

Starrs and Bartlett [io| have observed a 2/3 reduction 
in single-particle moduli for frequencies above 10 rad/s 
in polystyrene solutions. Using the solvent viscosity of 
decalin (2.6 mPa— s), a trap spring constant of « 1.66 x 
10-6 N/m, and Fig. 6 and 7 of Ref. [l^, r]uj/\ti\ can be es- 
timated to be 0.3 - 0.6 for frequencies between 10 rad/s 
and lO'' rad/s. To obtain 1/-^ ~ l/a for the semidi- 
lute polystyrene solution (concentration c = 1.1c*), we 
use a bead size a = 0.64 ^m [l3|, and estimate the cor- 
relation length as I = Rg{c*/c) w 0.06 /Ltm [1^, where 
Rg « 100 nm |13] is the radius of gyration. We obtain 
the estimate l/^/y ~ 10^^. From these estimates we 
conclude that uj may be slightly greater than lUs , but not 
by orders of magnitude. Therefore the effects of sliding 
may be important. 

On the other hand, Buchanan et. al. [4] have reported 
that the generalized Stokes Einstein relation is accurate 
in wormlike micellar solutions. Using the solvent viscos- 
ity of water, and the moduli from Fig. 1 of Ref. [j|, 
it can be estimated that riuj/\iJL\ is between 10"'^ and 10° 
for frequencies where the one particle microrheology data 
is reported, and in most cases ■qLu/\iM\ > 10"^. Using a 
mesh lengthscale of 10 nm and particle radius of 1 /im, 
1/^/7 Ki l/a « 10~^. Our calculation would then lead to 
the conclusion that for most of the data reported, and 
certainly for frequencies ly > 100 Hz [for which r]Oj/\iJ,\ is 
at least an order of magnitude greater than 10~^], the ef- 
fects of sliding should not be important. In this system, 
discrepancies from the generalized Stokes-Einstein rela- 
tion due to sliding effects can only appear for frequencies 
below « 1 Hz. 

In the examples provided above, we have focused on es- 
timating LUs, rather than LUf. However, if w/ < lUs, sliding 
effects will not appear until to < ivf. This dependence on 
friction may provide a way to test the predictions of our 
analysis by varying the surface chemistry of the probe 
particles @. 

While in most of this discussion we have focused on 
probe particles larger than the mesh size (7 « a'^/l'^ > 1), 
our results also apply when the probe particle is smaller 
than the mesh size (7 < 1). In this case, uj^ ~ ImI/('77) 
can be large and may be more easily accessible to ex- 



periments. In addition, when probe particles are much 
smaller than mesh size, they may be more likely to not 
directly interact with the network, leading to small val- 
ues of the friction coefficient and large values for LUf. Our 
solution provides a framework to analyze experiments in 
this regime of small probe particle size. 

Finally, an additional observation from the form of our 
results is that there is the possibility for features of the 
frequency dependence of moduli, including power laws, to 
be obscured or distorted in the vicinity of the crossover 
frequencies uJs, oJc, and ojf. This suggests that caution 
may be in order near these crossover frequencies. 

In the literature there has been much discussion of how 
probe-material interactions can affect one-particle mi- 
crorheological measurements. For example, near a probe 
particle the polymer network can be depleted by steric 
hindrance. One effect of this depletion zone is to change 
local rhcological properties, but another may be to in- 
crease the likelihood that the bead slides relative to the 
polymer network [lO| . Our calculation takes into account 
the second effect but not the first. Quantitative compari- 
son of Goff to experiments may require treatment of both 
effects, as well as the effects from adhesion of probe sur- 
face to polymer networks 3|. 
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APPENDIX A: DETAILS OF SOLUTION 

The solution to the equation for the stream function 
X can be obtained by using the axisymmetry of the so- 
lution, so that X is independent of (j). In addition, due 
to the boundary conditions at infinity, x raust be pro- 
portional to cos 0. Using these conditions in Eq. [^11 we 
find that the radial dependence of x satisfies an ordinary 
differential equation. The solution to x which results in 
finite displacements at infinity is 



X = Rc 



^ + i?r + J(l + ^)e-'=-'- 



(Al) 
where fc-^ = ^—i^p, in which the root with positive real 
part is chosen. 

To obtain Eq. [53] we need to find the matrix M which 
diagonalizes the equations for the stream function, 
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In the above time derivatives have been evaluated with 
respect to the oscillatory function exp(— ii) which multi- 
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plies all parts of the solutions. Explicitly, M is 
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1 
-1 
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(A3) 



and the eigenvalues are Aq — 0, A/j 
The solutions to a and (3 (Eq. 



(l-ii?)7. 
are found by a 

I 



similar procedure as the solution to x- ct and /? must 
be proportional to sin^ 9 to satisfy the boundary condi- 
tions. The radial dependence obeys an ordinary differ- 
ential equation, and the solutions which result in finite 
velocities and displacements at infinity are 



Re <^ Ar 



B 



Cr sin^ 



13 = Rej('i?r2 + :^+F(i + y^)e- 



(A4) 

(A5) 



Again, the square roots denote the root with positive real 
part. Note that it would not be difficult to generalize this 
method of solution to incorporate the effects of inertial 
terms. 

The boundary conditions Eq. [17] at r = a yield four 
(linear) equations for the coefficients A — G, since there 
both the velocity and displacement fields each have an r 
and 9 component. 

Additional conditions on the coefficients A — G arise 
from the boundary conditions for the displacement and 
velocity field at r — oo. Athough each have an r and a 9 
component, the r and 9 component give identical equa- 
tions for the coefficients A—G, leaving only 2 independent 
conditions out of the four: 



,r-i, 



2Affi^A 



2Aff2 £> + H 



(A6) 



J 



2M21 A + 2M22D 



-ib. 



(A7) 



In these equations the matrix elements of M ^ appear. 



Eq, 
exp(- 



M has 



terms 
and r 



proportional to exp(— A:^?'), 



^— A^r), r°, and r~^. The parts proportional to 
the exponentials are automatically satisfied by the gen- 
eral solutions in Eqs. lAll and lXsl The terms proportional 
to r*^ and r~^ give nontrivial relations for the coefficients 
A—G. The r and 9 component give identical equations 
for the coefficients, again leaving two independent condi- 
tions out of the four: 



i(2Mfi^A + 2M^^D + H) + {2M:^^ A + 2M:^^ D) 
i{2M^^B + 2M^^E - 2G) + {2M^i^B + 2M^2E) 

I 
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However, note that the first of these (from r^) is not in- 
dependent from the equations arising from the boundary 
conditions at r = 00. 



APPENDIX B: PROPERTIES OF THE 

RESPONSE FORCE WITH COMPRESSIBLE 

NETWORK 



Together, these produce only seven linearly indepen- 
dent equations for the nine coefficients A- J, so solving 
them leaves two coefficients undetermined. However, 
these undetermined coefficients only appear when the 
stream and potential functions are considered indepen- 
dently of one another; in the physically meaningful com- 
binations expressing the velocity and displacement field 
no undetermined coefficients remain. Since the solutions 
are not illuminating, we do not write them down. 



In this appendix we describe the properties of the re- 
sponse force when effects from compressibility of the net- 
work are included. First we isolate the effects of com- 
pressibility by examining the case of no-slip boundary 
conditions between the sphere and the network. Then 
we include effects from sliding of the sphere past the net- 
work. 



1. No-slip limit with compressible network 

In this section we address the hmit of no-sUp be- 
tween the sphere and the polymer, ^ — > cx) when A is 
finite. Examining the form of the solutions, one can see 
that the stream function x corresponding to the com- 
pressive motion of the network involves an exponential 
exp(— ry^P7/a) with length scale Re[a/y^p7]. The com- 
pressive motion dies off at longer scales due to the cou- 
pling to the incompressible fluid. Associated with these 
compressive effects is a crossover frequency Wc set by the 
condition |p|7 = 1, or ujc ~ |2/i -I- A|/(r/7). The expres- 
sion for ujc in the introduction is written in terms of vari- 
ables which can be easily controlled experimentally using 
the approximation 7 « a^//^, and hence applies when 
r « Tj/P. For frequencies less than Wc, compressive ef- 
fects of the network are important, while for frequencies 
above Wc, compressive effects of the network are not im- 
portant. This crossover is the same as that identified in 
previous studies [2, Q • 

For frequencies w 3> Wc, the real part of the modulus 



J 



G"ff(w < tJc) = w 1 + 



12 



GJjff tends to /ii , in accord with the expectations of Eq. [T] 
For frequencies lj -^ lOc, the real modulus is diminished, 
and for very low frequencies tends to the value 
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(Bl) 
Eq. IB II ranges from 80-100% of the high frequency limit 
/ii. For an incompressible network (A ^ 00) there is 
no reduction, while maximal reduction (20%) occurs for 
A = and /i2 ^ A*i. For /ii « ^2 ~ A (reasonable values 
for an actin network), the reduction is about 14%. In 
Fig. [5] the behavior of G^g is plotted. 

Now we turn to the imaginary part of the effective 
modulus G"ff. As in the incompressible case it is useful 
to analyze contribution with /i2 = separately from the 
polymer network contribution. At high frequencies, for 
/i2 = 0, G"ff — —ojr], in agreement with Eq. [T] For fre- 
quencies below ujc, G"ff crosses over to behavior in which 
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Note that as u) ^ 0, G'^g ~ uj^^"^. This behavior is plot- 
ted in Fig. [5^, which shows the quantity G"g/{r]uj) (an 
apparent viscosity), so that the high frequency regime 
appears as a constant, while the low frequency regime 
diverges as w"^". The low frequency divergence of 
the apparent viscosity results from the fact that at low 
frequencies, the compressive length scale Re[a/^/p7] di- 
verges, and so there is a larger and larger volume with 
appreciable relative motion between the solvent and net- 
work. The dissipation at low frequencies is increased 
by the friction between compressive modes of the poly- 
mer network and the incompressible solvent. To show 
that the crossover frequency is Uc, in one of the insets 
we plot G"g rescaled by its low frequency magnitude, 

{G':^ - w)/[G;'ff(^ « ^e)u=.. - V^] + 1- 



limiting value 



Now consider the case when the imaginary part of the 
network modulus {112) is nonzero. In addition to the con- 
tribution to G"ff described in the previous paragraph, 
there are also contributions from the network modulus 
/X2- For high frequencies (w > 10^), this additional con- 
tribution is simply 112 , as expected from Eq. [T] For low 
frequencies lu <ti ujc this contribution is diminished to the 
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As in the real part, the low frequency value is between 
80% (for A = 0, ^2 ^ fJ-i), and 100% (incompressible, 
A — 00) of the high frequency value. For /ii ~ /i2 ~ A, at 
low frequencies the low frequency additional contribution 
is about 81% of the high frequency additional contribu- 
tion. 

For a constant /i2, at high frequencies the contribution 
proportional to rjuj dominates, while at low frequencies 
the contribution from //2 dominates. This yields G'^g 
which behave as shown in Fig. [^b and c. Note that the 
low frequency enhancement discussed in the case /i2 = 
is not readily apparent since the contribution from /i2 
dominates at low frequencies. 



2. Frictionless limit with compressible netw^ork 

Now we turn to the opposite limit of no friction be- 
tween the polymer and sphere (^ = 0). As in the in- 
compressible case, we expect that for frequencies below 
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FIG. 8: (color online) Logarithmic plots of dimensionless G^ff 
as a function of dimensionless frequency for no-slip boundary 
conditions between sphere and polymer network. Eq. [T]would 
predict Gcfj = fj,\ . The shape of the plots indicates a crossover 
frequency ujc ~ |2/i + \\/(r)'y). For the effective modulus 
crosses over to new behavior in which compressional modes 
of the network are important. At low frequencies G^uj < uoc 
effects from compression cause G'^^ to underestimate /ii by 
up to 20%. b) Gpff for a range of 7 and \ = yL\ — fi2- 



(jjs the far-field solution takes the character of a solution 
with frictionless sliding boundary conditions. For fre- 
quencies uj > uig, the modulus GcS obeys Eq. [2 so that 
Gpff = Hi and G'^g = {r]Lu + H2)- For frequencies less than 
LUs, now that there is also compressibility, there are two 
cases to consider. In the first case, lUc < Ws, while in the 
second case ujc > i^s- 

First, consider the case where ouc < (^s- For Wc <C w ^ 
LUg, the real modulus G'^f^ = 2^i/3, 66% of the value in 
Eq. [TJ For lj <C Wc, compressive effects come into play 
and further diminish G'^g, which tends to 
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The low frequency limit of G'^g is between 57% (com- 
pressible limit with A ^ 0, /i2 ^ /^i) and 66% (incom- 
pressible limit, A ^ 00) of the high frequency value. For 
/i]^ « /Z2 ~ A, the low frequency limit of G'^g is about 
60% of the high frequency value ^i . 




rfuj/^ii 



b) 

ileff 
Ml 

10.0 
8.0 



1.0 
0.8 









// 


A=/li 7=1 ^2=/^! 

A=/ii 7=1 ^2=10/11 
A=8/ii 7=1 A'2=/'i 


A=^i 7=10 ^2=^*1 
\=Hi 7=10 /i2=10/Ji 
A=8/J.i 7=10 ^2=Aii 
A=0 7=1 A'2=10/ii 
A=0 7=10 /J.2=10^ti 


.^.^^^^^^ / 


/ 




»-"*^^ 


A \ 1 1 1 1 1 1 ll \ 1 1 1 1 1 1 ll \ 1 1 1 1 1 III \ 1 1 1 1 1 III \ — ^ 



10" 



10" 



10" 



iCr/cj/^iilo' 




10 ^^ r]uo/ Hi 



FIG. 9: (color online) Logarithmic plots of the imaginary part 
dimensionless G"fj as a function of dimensionless frequency 
for no-slip boundary conditions between sphere and polymer 
network, a) G"g/{rju}) as a function of frequency for 112 = 0. 
This contribution to G"ff from solvent viscosity and network- 
solvent drag dominates for small /i2 or high frequencies. At 



low frequencies [lj < cOc), G"fj 



,1/2 



(top inset). Lower 



inset: G^ff rescaled by its low-frequency magnitude (see text) 
to demonstrate that the crossover frequency is Uc- b) G"ff as a 
function of frequency. Eq.[T] would predict G"ff — fi2+riuj. For 
uj < uic, effects from compression cause G"ff to underestimatee 
fi2 by up to 20%. The crossover behavior is the same as in 
the real part except that for ui < H2/V the solvent viscosity 
(riuj) dominates, c) G"fj for a range of 7 and \ = Hi = H2- 
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In practice the condition that ujc <^ io <S^ Ws may not 
be met, either because LUg < lUc, or because LUf, is not 
sufficiently smaller than LUg. In this case G'^g appears 
to immediately tend to the behavior of Eq. IB4I when 
Lu < ujs. We plot representative behavior of G'^q in the 
frictionless limit in Fig. [TOl 



FIG. 10: (color online) Logarithmic plots of dimensionless 
Ggff as a function of dimensionless frequency for frictionless 
boundary conditions between sphere and polymer network. 
Eq. [l] would predict Ggg = ni. a) For frequencies less than 
both uJs and u>c, effects from sliding and compression cause 
G'gff to underestimate ni by up to 43%. b) G'^g for a range 
of 7 and X = fj.i = fj.2. li loc <^ uJs, ioT luc ^ io ^ lus, effects 
from sliding can cause G^g to underestimate /ii by 33%, and 
for w <C Wc <C Ws, compressional effects cause an additional 
underestimation by up to 43% of fii . 



Now we turn to the imaginary modulus G'^g. Again, 
first look at the contribution arising from the solvent vis- 
cosity and drag between the polymer network and sol- 
vent, by setting /i2 = 0. When w ^ Ws, the imagi- 
nary modulus G'"fj = —rjuj. For frequencies lu < lu^ 
crosses over to behavior of 
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The behavior of G"jf with /i2 = is shown in Fig. [TTk . 
which shows the apparent viscosity G"g/(?7Lij). In one of 
the insets we plot G"jj rescaled by its low frequency value, 
(G"ff-?7w)/[G"fj(u; < '^s)\uj=i^,-r\io\^\, to show that the 
crossover value is tOg- Note that for certain parameters, as 
Lo decreases, one can see the effective viscosity approach 
a plateau, which has the same origin as the plateau in the 
incompressible case described in the main text. At even 
lower frequencies, Eq. IBSl is dominated by w^' ^ behavior. 



This has the same origins in the compressibility of the 
network as discussed in the previous section dealing with 
no-slip boundary conditions. 

If /^2 7^ there is an additional contribution to the 
imaginary modulus. For tj 3> tj^, this contribution is /i2, 
the result of Eq.[T] When cjc <C w <C w^, this contribution 
is 2/i2/3. Finally, when lo < u^^ this contribution tends 
to the limiting value 
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FIG. 11: (color online) Logarithmic plots of dimensionless 
G"ff as a function of dimensionless frequency for frictionless 
boundary conditions between sphere and polymer network, a) 
G'cft/iv^) ^s a function of frequency, for ^2 ~ 0. This contri- 
bution to G"ff from solvent viscosity and network-solvent drag 
dominates for small fi2 or high frequencies. At the lowest fre- 
quencies G"ff ~ tj^" (left inset). Right inset: G"g rescaled by 
its low frequency magnitude (see text) to demonstrate that 
the crossover frequency is uja . b) Eq. [T] would predict that 
G"ff = jj,2 + rjuj. For frequencies less than both lJs and ujc, 
sliding and compressional effects cause G"ff to underestimate 
IJb2 by up to 43%. The crossover behavior is the same as in 
the real part except that for u) > ^J,2/r] the solvent viscosity 
(r/w) dominates, c) G"g for a range of 7 and X — ^\_ — 1^2- 
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This limiting value is smaller than the contribution from 
/i2 at high frequencies. It is between 57% (compressible 
limit, A = 0, /i2 ^ Hi), and 66% (incompressible limit, 
A — > oo) of the high frequency contribution correspond- 
ing to Eq. [T] For /ii « /i2 ~ A the limiting value is 58% 
of the high-frequency contribution. 

If LOc > ujs, for u; < Us the additional contribution to 
Ggg immediately tends to the value in Eq. 



3. Intermediate friction 

For any finite value of ^, there will be a crossover from 
behavior similar to ^ = at low frequencies to behavior 



similar to f ^ cx) at high frequencies. 



The crossover occurs at w/ — | =^ | . Typically, as plot- 
ted in Fig. [T^ because the no-slip hmit of G'^g is closer 
to Hi, G^jf appears to crossover from a low frequency 
limit of 2/ii/3 (or less, if compressional effects are impor- 
tant) to a value closer to the high frequency limit /ii at 
uj = imTi{ujf,ujs}. Similar behavior is seen in the imagi- 
nary modulus Ggg, shown in Fig. [121 
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FIG. 12: (color online) Logarithmic plots of dimensionless 
s-) Geff, b) G"fj as a function of dimensionless frequency 
for boundary conditions with intermediate friction between 
sphere and polymer network. In these plots \ = fj,\ = ^2 
and 7 = 3. The behavior interpolates between that of no- 
slip boundary conditions at high frequencies and frictionless 
boundary conditions at low frequencies, with crossover fre- 
quency u)f. For frequencies less than ujs, uif, and u)c, the 
moduli can be underestimated by up to 43%. 



